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Abstract. We introduce a Hartree-Fock-Bogoliubov mean-field approach to treat the problem of proton 
emission from a deformed nucleus. By substituting a rigid rotor in a particle-rotor-model with a mean-field 
we obtain a better description of experimental data in 141 Ho. The approach also elucidates the softening 
of kinematic coupling between particle and collective rotation, the Coriolis attenuation problem. 
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Proton emission is a weak single-particle (s.p.) process 
with widths about 20 orders of magnitude smaller than the 
usual MeV scale of other nuclear interactions. This makes 
observation of proton radioactivity an ideal and power- 
ful tool for non-invasive probing of the single-proton in- 
medium dynamics. Recent studies have already explored 
numerous nuclear mean-field properties of proton emitters 
including deformations, vibrations rotations 0, pairing 
and other many-body correlations 

In this work, using proton emission from deformed nu- 
clei, we concentrate on an old problem known as Cori- 
olis attenuation problem in the particle-rotor model 
(PRM). Recent studies of proton decay 2 4 highlight the 
same lack of kinematic coupling between the particle and 
the deformed rotor as was inferred decades ago from ob- 
servations of the energy spectra of odd- A nuclei The 
second purpose of this work is to gain an understanding 
of and to develop a better theoretical technique to de- 
scribe particle motion in the deformed mean-held. Here 
the notion of a core as a rigid rotor is inadequate and, as 
emphasized in numerous works (Sfl7||5], the residual two- 
body interaction and collective modes are important parts 
of the dynamics. 

We consider an axially-symmetric deformed proton emit- 
ter and assume that the total Hamiltonian is composed of 
a collective H co n = R 2 L /2 J C and intrinsic parts 
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Here R denotes the rotor angular momentum, involving 
only the part perpendicular (_L) to the symmetry axis, and 
o/q and a n stand for s.p. creation and annihilation opera- 
tors of state \f2) in the deformed body- fixed mean-held 
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potential. Nuclear pairing involves body-fixed time(7?.)- 
conjugate s.p. states \f2) and \S7) and describes the resid- 
ual two-body interaction. In contrast to the usual PRM 
this model assumes some odd number of valence parti- 
cles. In the limit where the valence space covers the entire 
nucleus the collective rotor variables become redundant. 

Kinematic coupling between the intrinsic system and 
collective rotor occurs due to conservation of total angular 
momentum I = R + j, where j is the angular momentum 
of the valence particles. Components of this operator can 
be expressed in the a intrinsic body-fixed basis as 
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similarly for = The coefficients jnsi 1 — (^b+l^') 
are obtained using expansion of states \ fi) in spherical 
basis. Excluding a trivial rotational part from the total 
Hamiltonian H = I 2 /(2£) + H' we obtain 
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which is to be solved via many-body techniques using basis 
states formed as products of Wigner £>{ /x (cj)-functions 
of collective angles lo, and any complete set of many-body 
intrinsic states such as Slater determinants. 

Here we implement a Hartree-Fock-Bogoliubov (HFB) 
approach that allows one to determine a s.p. mean-held, 
which is a combination of the rotor degrees of freedom and 
even-particle valence system, and absorbs in the best way 
kinematic couplings and residual nucleon-nucleon correla- 
tions. By making a Bogoliubov transformation to quasi- 

particles — J2n ( u h a Si + v h a Sij anc ^ with the require- 
ment that the elementary quasiparticle excitations are sta- 
tionary we obtain the usual HFB equations 
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Here in full analogy to PRM the diagonal part of the s.p. 
potential is given by the usual s.p. energy corrected with 
the recoil term and decompiling factor AE jS] 



enn = e n + ^{{n\i 2 \Q)-2Q + 5 Shl/2 AE) . (5) 



The off-diagonal term in Eq. (jJJ violates deformation align- 
ment, the iGsymmetry, which manifests itself through 
non- vanishing average mean-field expectations = 
(j) while (js) = 0. This average mean-field value enters the 
off-diagonal s.p. potential 
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and is to be determined in a self-consistent solution 

(?) = 2 3n+i,{2 v l Q+l v l n . 
i, n>o 
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This is analogous to non-conservation of particle number 
N, a common situation in the HFB approach. Particle 
number is restored on average via the introduction of a 
chemical potential H' — > H' — fiN, so that the pairing 
gap and chemical potential in Eq. are self-consistently 
determined 
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The term (j) in Eq. JBJ is due to HFB linearization 
of the recoil operator j 2 ~ + J-)/2 + f2 2 which, 

besides acting on an odd particle, also perturbs an even- 
particle mean-field, thus producing a suppression of the 
Coriolis mixing. The Coriolis interaction takes the form 
— (I — similar to the Routhian in the Cranking 

Model Uj|, and is suppressed. This is in contrast with the 
PRM, where by definition the rotor is rigid and (j) = 0. 

The quantity £= (l- - is the av- 

erage suppression factor; for the case of 141 Ho (see below), 
it is shown as a function of pairing gap in the figure. The 
idea to phenomenologically substitute the spin of the rotor 
R = (I— for the operator I in order to explain Coriolis 
attenuation was suggested in 0, and contributions from 
the j 2 operator in the mean-field approach are discussed 
in Other contributions coming from non-rigidity of 

the core are also considered j3Hj- 

We apply this approach to the proton emitter 141 Ho 
where partial decay widths Jo for decay to the + ground 
state and F2 to the 2 + first excited state in 140 Dy are 
known from experiment. The spectrum of 140 Dy is used 
to determine deformation and moment of inertia. The va- 
lence space is limited to a negative parity subspace coming 
from spherical /in/2 orbital, but particle depletion due to 
pair excitation onto positive parity states is included. The 
decay amplitudes computed using appropriate deformed 
Woods-Saxon potential and expressed via normalization 
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Experiment 


10.9 
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of the wave function [2] Apjk) = <j)?Ar) / 'Gy(fcr) 
where Gij is the irregular Coulomb function. The decay 
width is given by 0] r 



k 2(2^+1) 



21+1 Z^f2>o ^jK,m a n Ji ij 
where C is a Clebsch-Gordan coefficient and the Uq fac- 
tors come from the solution of Eq. 0} . The results of this 
calculation, labeled as RHFB, are compared with PRM 
and experiment in the Table. The problem with Corio- 
lis attenuation in PRM is transparent in this case; e.g., 
for Jo (first column), the adiabatic limit (C — > 00) over- 
estimates experiment, while introduction of Coriolis mix- 
ing even softened by pairing correlations extremely over- 
reduces To- The HFB calculation shown in the Table, al- 
though limited to a very small valence space, already leads 
to a substantial improvement. 
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